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Abstract 

| Subcritical catalytic branching random walk on (f-dimensional lattice is studied. New the- 

orems concerning the asymptotic behavior of distributions of local particles numbers are es- 
tablished. To prove the results different approaches are used including the connection between 
fractional moments of random variables and fractional derivatives of their Laplace transforms. 
■ In the previous papers on this subject only supercritical and critical regimes were investigated 

Q . assuming finiteness of the first moment of offspring number and finiteness of the variance of 

offspring number, respectively. In the present paper for the offspring number in subcritical 
regime finiteness of the moment of order 1 + S is required where 5 is some positive number. 
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spring number, infinite variance of offspring number, local particles numbers, limit theorems, 
fractional moments, fractional derivatives. 



43 



2010 AMS classification: 60F05. 



1 Introduction and main results 

oo 

^ The paper is devoted to investigation of catalytic branching random walk (CBRW) on integer lattice 
Z d , d G N. This modification of branching random walk (BRW) with a single source of branching 

i— h was proposed by V.A. Vatutin, V.A. Topchij and E.B. Yarovaya in [1] and it comprises symmetric 

^ - BRW studied earlier (see, e.g., [2]). 

Recall description of CBRW on Z d . Let at the initial time t — there be a single particle on the 

.£h. lattice located at point x G Z d . If x ^ then the particle performs the continuous time random 

^ walk until the first hitting of the origin. We assume that the random walk outside the origin is 
c3 specified by infinitesimal matrix A = (a(u, v)) u ve zd and is symmetric, homogeneous, having a finite 
variance of jumps and irreducible (i.e. the particle passes from an arbitrary point u G Z d to any 
point v G Z d within finite time with positive probability). It means that 



a(u, v) = a(v, u), a(u, v) = a(0, v — u) := a(v — u), u, v G Z d , 
?d a(v) = where a(0) < and a(v) > for v^O, d ||v|| 2 a(v) < oo 



(here || • || denotes an arbitrary norm in M. d ). If x = or the particle has just hit the origin then it 
spends there random time distributed according to exponential law with parameter 1. Afterwards 
it either dies with probability a, producing just before the death a random offspring number £ or 
leaves the origin with probability 1 — a, so that the intensity of transition from the origin to point 
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v / equals —(1 — a)a(v)/a(0). At the origin the branching of particle is specified by probability 
generating function 

f(s):=EsZ = Y7 JkS k , se[0,l}- 

*— »fc=0 

At the birth moment the newborn particles are located at the origin. They evolve in accordance 
with the scheme described above independently of each other as well as of their parents history. 

The natural objects of study in CBRW are total and local particles numbers. Denote by fj,(t) 
the number of particles existing on the lattice Z d at time t > 0. In a similar way we define local 
numbers fj,(t; y) as quantities of particles located at separate points y G Z rf at time t. 

It was established in [3] that exponential growth (as t — > oo) of both total and local particles 
numbers in CBRW on Z d holds if and only if E£ > a). Here hd is probability of the event 

that a particle which has left the origin will never come back. Thus, the value E£ = l + — a) 

is critical and, similarly to many types of branching processes (see, e.g., [1]), CBRW is classified 
as supercritical, critical or subcritical if the mean offspring number E£ is greater, equal or less 
than 1 + hdOT 1 ^ — a), respectively. In view of the random walk properties such as recurrence or 
transience, h± = /t2 = 0, whereas < hd < 1 for d > 3. Recall that E£ = f'(l) and for classic 
Galton- Watson branching processes the critical value of f'(l) is equal to 1. 

Critical and subcritical CBRW on Z d are of special interest since in these cases there arise diverse 
kinds of limit (in time) behavior of the total and local particles numbers depending on dimension d. 
For instance, for d — 1 and d = 2 probability P x (/x(£) > 0) of non-extinction of the population tends 
to zero as time grows (index x G Z d denotes the starting point of CBRW), whereas for d > 3 this 
probability has a positive limit. Such effect is due to existence, on lattice Z d for d > 3, of particles 
which with positive probability are ever alive and never hit the source of branching. The total 
particles number in critical CBRW on integer line Z, i.e. for d = 1, was studied in the fundamental 
paper pQ. In papers [5| - [8 J the investigation was continued for critical and subcritical CBRW on 
7L d for any d G N. 

The analysis of local particles numbers is much more hard. For critical CBRW on Z d the 
limit distributions of local particles numbers were studied in the series of papers by V.A. Vatutin, 
V.A. Topchij, Y. Hu, E.B. Yarovaya and the author (see, e.g., [H] - [12] )• It is necessary to note 
that in all these papers the finiteness of second moment E£ 2 of offspring number was assumed. In 
the present paper the limit distributions of local particles numbers in subcritical CBRW on Z d are 
studied for the first time and the corresponding conditions on the moments of offspring number 
£ are less restrictive than those in the critical case. Namely, we find the asymptotic behavior 
of the mean local numbers m(i;x, y) := E x /z(i; y) for fixed x, y G Z d and t — > oo. Assuming 
that E£ 1+5 < oo for some 5 G (0,1], the similar problem is solved for non-extinction probability 
q(t; x, y) := P x (/x(t; y) > 0) of local particles numbers. Moreover, under the same restriction on the 
moment E£ 1+<5 the conditional limit theorems are proved for y) as t — > oo. 

To formulate the main results let us introduce additional notation. Let 

POO 

g(s,t;x,y) := 1- E xS ^ ;y ) and J(s;y) := / $(g(s, t; 0, y)) dt 

Jo 

where $(s) := ot(f(l — s) — 1 + f'(l)s) and s G [0,1], t > 0, x, y G Z d . Consider transition 
probabilities p(t;x, y), t > 0, x, y G Z d , of the random walk generated by matrix A. According to 
[H] and Theorem 2.1.1 in [T3], for fixed x and y, as t — > oo, the following asymptotic relations hold 
true _ 

Pi^y)-^ P'teM)--^^' P(*;0,0)- P (t;x,y)~ 7 ^ 1 x) (1) 
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where ld := ((2vr) rf |det 0^(0) |) 1/2 , 0(0) : = £ zeZd o(0, z) cos(z,0), G 
& 



:= 



, ^(z):=-— -/ (v,z)V^°) v ' v >dv, ZGZ" 

0/ ij6{l,...,d} ^7TJ 



and (•, •) denotes the inner product in R d . Set G\(x, y) := J °° e~ xt p(t; x, y) dt, A > 0, x, y G Z d , i.e. 
G^(x, y) is the Laplace transform of the transition probability x, y). By virtue of ([T]) for d > 3 
the Green's function Go( x ; y) is finite for all x, y G Z d , which means the transience of the random 
walk. However, for d — 1 or d — 2, we deal with recurrent random walk since liniA^o+ Ga(x, y) = oo. 
In view of the same formula (00) the function liniA^o+ (G\(0, 0) — Ga(x, y)) is finite for all d G N 
and x, y G Z d . So we may introduce the function 



I (l-a)a- 1 -Pj °°(p(t;0,0)-p(t;0,z))dt if z ^ 0, 
MzJ-i j if z = 



^(0,y):= 7r ^p 2 (y), C 2 (x,y):= — -p 2 (x)p 2 (y), x ^ 0. 



where for the sake of convenience we set a := — a(0) and /3 := a(f'(l) — 1). As shown in [9], 
/id = (aGo(0, 0)) _1 . It implies that in subcritical regime (3 < (1 — a)(aGo(0, 0))" 1 and, therefore, 
Pd(-) is a strictly positive function for all d G N. Note also that in accordance with formula (2.1.15) 
in [T3] inequality 0, 0) > p(i;x, y), t > 0, is valid and the function p(-;x, y) is symmetric and 
homogeneous in x, y G Z d . 

For d = 1 and x, y G Z, we define 

1 - a 1 

Ci(0,y :=- ^Pi(y), Ci(x,y :=- Pl x Pl (y + 7l x) + 7l (y - 7l (y - x , x ^ 0. 

za 7 i7rp z 2 7 l7Tp z 

In a similar way, for d = 2 and x, y G Z 2 set 

^|p 2 (y), C 2 (x, y): =-i 

At last, for d > 3, specify functions Cd(x, y), x, y G Z d , by way of 

Cd(0, y) := 7: — a )^Z d n 7^ Pd{y), C d (x,y) := p d (x)p d (y), x ^ 0. 

(1 — a — apG (0,0)J' ! (1 — a — apGo(0,0)) 

Theorem 1 Let E£ < l+/i d af~ :L (l — a). T/ien, ast — > oo, for each x, y G Z d ; the following relations 
are true 

w(t;x,y) -p — , a = 1, 

/. x C2(x,y) 

m{t;x,y) 2 , d = 2, 

fin t 

C d (x,y) 

wK*;x,y) — , rf>3, 

and f/je functions C d (-, ■), d G N, introduced above, are strictly positive. 

The statement of Theorem [T] generalizes the corresponding results of Chapter 5 in [13] concerning 
the asymptotic behavior of the first moments of local particles numbers in subcritical symmetric 
BRW on Z d . Here, in contrast to [TBj, we use the approaches such as representation of complex- 
valued measures in terms of Banach algebras (see [9]) and Tauberian theorems for derivatives of 
Laplace transform (see Section 7.3 in [14|). 
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Theorem 2 If E£ < 1 + h^a l (l — a) and there exists E£ 1+<5 for some 5 G (0, 1] then, for fixed 
x, y G Z d ; as t — > oo, one has 

( , v Ci(x,y)-C 1 (x,O)J(0;y) 

( . v C 2 (x,0)(p 2 (y)-J(0;y)) 
9(*;x,y) , d = 2, 

, g d (x,0)(p d (y)- J(0;y)) 

?(*; x > y) ~ , d > 3, 

where Ci(x, y) — Ci(x, 0) J(0; y) > /or a// x, y G Z and J(0; y) < Pd(y) /or d > 2 and y G Z' 



To prove Theorem[2]we apply the Holder inequality combined with results on connection between 
fractional moments of random variables and fractional derivatives of their Laplace transforms (see, 
e.g., [15J). Theorem [3] can be considered as a corollary of Theorem [3 

Theorem 3 Let E£ < 1 + hdd^il — a) and E£ 1+<5 be finite for some 5 G (0, 1]. Then, for fixed 
x, y G Z d and each s G [0, 1], the following equalities are valid 

Urn E x (^\ Mi; y) > 0) = ^■(x.y)-C,(x,0)(J(0; y )-J( S ; y )) rf = 
1 Ci(x,y) - Ci(x,0)J(0;y) 



t— >oo 



iimE x ( s ^)|^;y)>o)= sMy) "'f ; ^f ;y)) . i>2. 

' - 1 Pd{y) - J(0;y) 



t—>-oo 



Remark. Comparing formulations of Theorems [TH3] and the results on local particles numbers 
in critical CBRW on Z 2 (see, e.g., [16] and [17]) the following conclusion suggests itself. Local 
particles numbers in critical CBRW demonstrate "subcritical" behavior in the case of random walk 
on integer lattice. Moreover, the scheme of proofs of Theorems [TH3] is easily carried over to the case 
of critical CBRW on Z 2 . Consequently, we may state that the results of [16] and [T7J, concerning the 
local numbers in critical CBRW on Z 2 , are valid under less restrictive conditions on the moments 
of offspring number. Namely, it is sufficient to require finiteness of E£ 1+<5 for some S G (0, 1] instead 
of condition E£ 2 < oo. 

Concluding the first part of the paper let us note the close relation between CBRW and su- 
perprocesses, namely, catalytic super-Brownian motion with a single point of catalysis (see, e.g., 
[18] and references therein). It is of interest that in view of [1] CBRW may be considered as a 
queueing system with a random number of independent servers. This gives an opportunity for wide 
applications of the established results. 



2 Proof of Theorem U 

Similarly to the proof of Theorem 1 in [T7J, we use backward and forward integral equations for 
the family of functions {m(-; x, y)} Xiye zd. These equations coincide with equations (8) and (9) in 
|17j . derived for the mean local particles numbers in critical CBRW on Z d , upon the replacement 
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of critical value /3 C = (1 — a)a 1 G 1 (0, 0) by value (3, that is, 

m(£; x, y) = p(t; x, y) + ^1 - ^ °_ - j ^ p(t - u; x, 0)m'(>; 0, y) cfcz 

+ t^— / p(* - w; x, 0)m(u; 0, y) du, (2) 

m(t; x, y) = p(t; x, y) + ^ - - - - lj ^ m(it; x, 0)p' (t - it; 0, y) dtt 

+ f3 m(u; x, 0)p(t — it; 0, y) cfat. (3) 
Jo 

Relations d2J) and (EJ), as well as equations (8) and (9) in [T7], are obtained by means of the 
variation of constants formula applied to backward and forward differential equations (5) and (6) 
in [T7] established in Banach space l oc (Ti d ). 

We also need the following auxiliary statement proved similarly to Lemma 3.3.5 in [13J and 
Lemma 1 in [17J. 

Lemma 1 For each y G 7L d the function m(t; y, y) does not increase in variable t. 

Now we turn directly to the proof of Theorem [TJ At first let us consider the case x = y = 0. 
Apply the Laplace transform to both sides of equality (j3J) and use the obtained relation to express 
function m(A) := J °° e" A *m(t; 0, 0) alt, A > 0. Then 

am = (4) 

K J 1 - ((1 - a)a- 1 - 1) J °° e- xt p'(t; 0, 0) dt - 0G X (O, 0) ' 1 ' 

Differentiate each side of the last relation in A. Consequently, taking into account the identity 
J °° e- xt p'(t; 0, 0) dt = AG A (0, 0) - 1 one has 



= (1 - cQq-^O, 0) + ((1 - cQa- 1 - 1)G|(0, 0) 
((1 - «)a~ 1 - ((1 - a)a~ l - 1) AG A (0, 0) - /3G A (0, 0)) 2 ' 



According to Tauberian Theorem 2 in §5 of Chapter XIII in [19] along with Corollary 43 in [H], 
relation (CQ) implies 



G A (0,0)~^, G A (0,0)~-^, d=l, 
G A (0,0)~ 72 ln^, G A (0,0)~-^, d = 2, 

as A — > 0+. Substituting these asymptotic equalities into (EJ) for d — 1 and d = 2, respectively, we 
find 

ffi'(X) ~ =, d — 1, 

2a 7lv /?/3 2 v / A 



1 — a 
a 72 /3 2 A ln"A : 



m'(A)~- — -^ttxt, d=2, 
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as A — > 0+. Applying Corollary 43 in [H] to the above relations we come to the assertion of 
Theorem [I] when x = y = and d = 1 or d = 2. 

For d > 3, we employ another approach, namely, the representation of complex- valued measures 
in terms of Banach algebras. In view of ([[]) two cumulative distribution functions, having the Laplace 
transforms G\(Q, 0) and f Q °° e~ xt p'(t; 0, 0) dt, respectively, possess tails equivalent to constants (the 
second constant being zero) multiplied by the same function t 1-d / 2 . Hence, due to Lemma 6 in [9] 
and formula (j3J) the following asymptotic equality holds true 

, 2 (1 — a) a^d 

m w; 0, 0) cm ~ — — — . . . , , n - , t — > oo. 

{d-2){\- a- a(5G Q {0,0)) 2 t d / 2 - v 



Whence by Lemma [T] and the classical results on differentiating of asymptotic formulae (see, e.g., 
[2D] . Chapter 7, Section 3) we get the assertion of Theorem [T] when x = y = and d > 3. 

Let us consider the case and y = 0. Integration by parts permits to rewrite the family of 
equations <^ as follows 

m(t;x,0)= Q p(t;x,0) + (l — ^ / m(t - u; 0, 0)p'(u; x, 0) 

1 - a V 1 _ a / Jo 

+ / pit — u; x, 0) m(u; 0, 0) du, (6) 

1 - " Jo 



— m(t; 0,0) = —^—p(t; 0,0) + fl — ^ / m(t - u; 0, 0)p'(u; 0, 0) du 

1 — a 1 — a \ 1 — ay Jq 

+ [ p{t-u;0,0)m(u;0,0)du. (7) 

1 - a Jo 



Subtracting equation (jSj) from (j7|) we come to 

m(t;0,0) -m(t;x,0) = — ^— (p(t; 0, 0) - p(t; x, 0)) 



1 — a 1 — a 

+ ( 1_ j - " ~ ) jf m(t - u; 0,0) (p'(u; 0,0) - p'(u;x,0)) du 

+ — ^- / m(t - it; 0, 0) (p(u; 0,0) -p(u;x,0)) d«. (8) 
1-ttJo 

Employing the results on differentiating of asymptotic formulae once again (see, e.g., |20j, Chapter 7, 
Section 3), as well as relation (P) and inequality p"(t; 0, 0) > p"(t; x, 0), t > 0, implied by (2.1.15) 
in [13], we find that 

pU0,Q)-^x,0)~- ( ^y , t^oo. (9) 

Therefore, on account of Lemma 5.1.2 ("lemma on convolutions") in [13] along with formula ([T]) 
and the proved part of Theorem [1] we deduce from ([8]) that 

m(t;x,0) ~ m(t; 0, 0) M - f (p(u;0,0) - p(u;x,0)) du 



(1 - a) t d l 2+1 



(^l + f3^ m(u;0,0)du S J , t ->■ oo. (10) 
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However by virtue of already established part of Theorem [T] we see that / °° m(u; 0, 0) du < oo for 
all d G N and, moreover, according to (@J for A = one gets 



/"OO 

/ m(u;0,0)du = -/T 1 if d = 1 or d = 2. 
Jo 



Hence we conclude that only the first summand in the right-hand side of fflO|) makes contribution 
to the asymptotic behavior of m(£;x, 0). So, the assertion of Theorem [T] for x / and y = is 
entailed by relation ( imp and the proved part of this theorem when x = y = 0. 
Let now xGZ d and y ^ 0. In view of d2J) one has 

m(t; x, 0) - m(t; x, y) = p(t; x, 0) - p(t; x, y) 



+ /3 m(t — u; x, 0) (p(u; 0, 0) — p(u; 0, y)) dtt. 



Then taking into account formulae (pj) and ([9]) along with Lemma 5.1.2 in [13] and the proved part 
of Theorem [1] we find that 

m(t;x,y) ~ m(t; x, 0) ~ P J (p(u; 0, 0) - p(u; 0, y)) duj 

+ 7^77 ^TTJT- / m( M ;x,0)rf M , ^ oo. (12) 



fd/2+l fd/2+1 

By the established part of Theorem [1] applied one again we come to inequality / °° m(u; x, 0) du < oo. 
Moreover, similarly to the verification of equality (1111) we check that J °° m(u; x, 0) du = —f3~ 1 for 
d = 1 or d = 2. Thus, the statement of Theorem [T] for x6Z d and y ^ is implied by relation (1121) 
and the proved part of Theorem [1] for x e Z d and y = 0. 

To complete the proof of Theorem [1] one has to make sure only that functions C<f(-, ■), <i € N, 
are strictly positive. It is easy except for the case d — 1 when x / and y / 0. Let us show 
that in this case C*i(x, y) > as well. For this purpose we turn to function H x (t), t > 0, which is 
a cumulative distribution function of time from leaving point x till the first hitting point in the 
framework of the random walk generated by matrix A. Clearly, 

p(*;x,y)- / p(t-u;O,y)dH x>o (u)>0, t>0. (13) 

Then by virtue of the evident identity p(t; x, 0) = f* p(t — u; 0, 0) dH x0 (u) one has 



P(*;x,y)-/ p(t-u;0,y)dH XtO (u) = p(t; x, y) - p(t; x, 0) 



+ / (p(t-u;0,0)-p(t-u;0,y)) dH Xt0 { 
Jo 



u). 



Since with the help of relation (jTj), Lemma 5.1.2 in [TJ] and Lemma 3 in [TT] one can find the 
asymptotic behavior of the right-hand side of the latter equality, we establish that, as t — > oo, 

7i(x) +7i(y) -7i(y-x) 



p(t; x, y) - / p(t - u; 0, y) dH Xi0 (u) 
Jo 



+ 



t 3/2 

[l - a- qpi(x))(l - a - api(y)) 
2a 2 7 1 7r / 3 2 t 3 / 2 
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Hence, it follows from ( ITSj) that 



Pi(x)pi(y) - / ~ / N ~ , 1 — a ( , x / \ I — a 

W + 7>(x) + 7:(y) - 7:(y - x) > ( ft (x) + p l( y) - — 

However the left-hand side of this inequality appears to be C*i(x, y), whereas the right-hand side is 
strictly positive since pi(z) > (1 — a)a~ x for z / 0. In its turn, the last inequality is satisfied due 
to definition of function pi(-) and negativeness of /3 in subcritical regime for d — 1. 
Therefore, Theorem [TJ is proved completely. □ 



3 Proofs of Theorems 2 and 3 



It is not difficult to verify (following the scheme in [T7]) that in subcritical CBRW on Z d , as well as 
in critical CBRW, for all x, y G Z d , s G [0, 1] and t > 0, the non-linear integral equations hold true 



q(s, t; x, y) = (1 - s)m(t; x, y) - / m(t - u; x, 0)$(q(s, u; 0, y)) du. 

Jo 



(14) 



So, the following upper estimate for g(s,t;x, y) ensues in view of non-negativeness of functions 
m(-; x, 0) and $(•) 

g(s,t;x,y) < (1 - s)m(t;x,y). (15) 

Lemma 2 // < 1 + /ida _1 (l — a) and E£ 1+5 < oo /or 5 G (0, 1] t/ien for some positive constants 
Ki and K 2 the inequalities are valid 

Hs) < K lS 1+s , SG[0,1], (16) 
E x p(t;y) 1+5 < K 2 m(t;x,y), t > t (x, y), x, y G Z d (17) 

with a certain non-negative function t (-,-). 

Proof. At first we consider the case < 5 < 1. Let us take advantage of the connection between 
fractional moments of random variables and fractional derivatives of their Laplace transforms. For 
the first time such results were obtained in [T5] where the traditional notion of Riemann-Liouville 
fractional derivative was used. However it is more convenient to involve the up-to-date counterpart 
of these results, namely, Lemma 2.1 in [21 J, which gives 



r(l — o) J av 2+d 

In view of finiteness of E£ 1+<5 the latter equality entails v~ 2 ~ 5 Q(v) dv < oo. Taking into account 
that $'(s) > 0, s G [0, 1], integration by parts implies relation f lT6|) for 5 G (0, 1). 

Turn to verification of f[T7|) when < 5 < 1. Applying Lemma 2.1 in [2TJ once again we come 
to the following equality 



(5(1 + (5) /-°° ?;m(t;x,y)-g(e-V;x,y) 
xM,yj ~T(1-5)J 
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Substituting formula (JUj) into the latter relation we get 



ExM^;y) + = m(t;x,y) ^^— ^ y ^ cfc 



5(1 + 5) r° i 



+ ^ ±1 I m(t-u-x,0)<5>(q(e- v ,u-0,y))dudv. (19) 



Obviously, in equality ( !T9|) the first integral converges. Let us estimate the double integral. Ac- 
cording to ( 1151) and inequality <&(/cs) < «$(s), k, s G [0,1], guaranteed by convexity property of 
function $(■), one has 



oo t t 



f 1 f _ /" /'$(l-e _1 ') 
/ ~~ 2+5 / m (^ ~ M ; x ) 0)$(g(e 0, y)) dwdi; < / m{t — u; x, 0)m(u; 0, y) <iu / ^ <it>. 

oo o o 

At the right-hand side of the above inequality the integral in variable u is equivalent (up to a 
constant factor) to the function m(t; x, y), as t — > oo, on account of Theorem [1] and Lemma 5.1.2 in 
[T3] . Furthermore, in the same inequality the integral in variable v converges by virtue of formula 
(|T8|) combined with finiteness of E£ 1+5 . Therefore, this arguing together with relation ( 1T9|) lead to 
(IT?]) with < 5 < 1. 

Now we consider the case 8 = 1. The identity f"(l) = E£(£— 1) and the Lemma conditions 
imply the existence of f"{l) and, consequently, the existence of $"(0). Moreover, since $(0) = and 
$'(0) = 0, inequality ( IPS"]) for 5 = 1 is proved. Let us take the second left derivatives at s = 1 for each 
side of equality ( TT4"1) . Using relations m(i;x, y) = — d a q(s,t',x.,y)\ a _ 1 and E x /i(t; y) (//(£; y) — 1) = 
- <9 S 2 S tf(s, t; x, y)| s=1 , we obtain 

E x /i(t;y)(/i(t;y)-l) = a/ ,/ (l) / m(t - w; x, 0) (m(u; 0, y)) 2 du. 

Jo 

In accordance with Theorem [1] and Lemma 5.1.2 in [13], in the latter equality the integral behaves 
as m(t] x, 0) up to a constant factor, as t — > oo. Hence, on account of Theorem [1] this entails the 
desired inequality ( 1171) when 8=1. Lemma [2] is proved completely. □ 

Let us turn to proving Theorem It is easily verified that due to formulae (ITBl) and ( fl~6l) along 
with Theorem [1] and the proof scheme of Lemma 4 in [16], one gets 

/ m(t - u; x, 0)$(g(s,w; 0,y)) du ~ m(t; x, 0) J(s; y), t ->■ oo. (20) 
io 

Observe that the function J(-; •), appearing in formulations of Theorems [2] and [31 is defined correctly 
in view of upper estimate (fT5|) for g(s,t;x, y). Let us find the lower estimate for g(t;x, y). By the 
Holder inequality one has 

E x Mt;y) = E x /i(t; y )!(//(*; y) > 0) < (E xM (t; y) 1+5 ) V(1+5) (E x I(^;y) > O)^)^^ 
where I(-B) denotes the indicator of set S. Rewrite the last inequality as follows 

(E x/ ,(t;y) 1+5 ) 1/5 
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Now employing (|17p we come to the desired lower estimate for q(t; x, y) 

g(t;x,y) >/^ 1/5 m(t;x,y), t>t (x,y), x,yGZ d (21) 
Combining formulae ( III"]) and ( )20l) . when s = 0, with estimate ( 12T1) we conclude that, as t — > oo, 

ttx (/■! x, y ) 

<?(£; x, y) ~ m(t; x, y) - m(t] x, 0) J(0; y) and J (0; y) < lim ' ' 

t->oo m{t; x, 0) 

for fixed x, y e Z d . Just these relations amount to validity of Theorem [2J □ 

Let us prove Theorem [3j Applying formulae tHM and (|20T) once again we have 

g(s,t;x,y) ~ (1 - s)m(t; x, y) - m(t; x, 0) J(s; y), t -> oo, x,y 6 Z d . 

Then with the help of the identity 

lim E x (s^)| ( t; y) > 0) = 1 - lim 

t^oo v 1 7 t^oo q(t; x, y) 

and Theorem [2] we obtain the assertion of Theorem [31 □ 
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